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Abstract —This paper presents the design and implementation 
of a novel SDR based massive MIMO testbed with up to 70 nodes 
built at Tennessee Technological University. The deployment can 
reach a 30 x 30 antenna MIMO scheme. With this testbed, we are 
able to measure the channel matrix and compute the achievable 
rate of the massive MIMO system using experimental data. The 
measured channel capacity is linearly increasing with the number 
of antennas of the base station. We also demonstrate the channel 
reciprocity including the circuits impact from the transmitter 
and receiver. We show that the Vandermonde channel model is 
more realistic to describe the massive MIMO architecture than 
the widely used Gaussian channel model, in terms of capacity. 
By adjusting the range for angle of arrival a and the base 
station antenna distance d during the simulation, we find out 
the Vandermonde model agrees with our measured capacity at a 
certain a for each selected d and the a is very close to that of the 
experiment deployment. It is the first time that the feasibility of 
Vandermonde channel model is demonstrated by the experiment 
for massive MIMO. 

Index Terms —Massive MIMO, Software Defined Radio, 
testbed, Channel Model 

I. INTRODUCTION 

Massive MIMO, as an candidate for one of the disruptive 
technologies of the next generation (5G) communications 
system, promises significant gains in wireless data rates and 
link reliability ID E). By using large numbers of antennas 
at the base station, massive MIMO helps to focus the radiated 
energy toward the intended direction while minimizing the 
intra and intercell interference. 

While promising, the research challenges of massive MIMO 
also emerge with the dramatically increased number of anten¬ 
nas at the base station. Obtaining channel state information 
is an essential operation for communications but the pilot 
contamination could occur due to the limited pilots for massive 
MIMO. 

Most of the current theoretical analysis on large MIMO 
assume a Gaussian Channel matrix m 0 However, Gaus¬ 
sian channel could be not accurate in real environment. The 
channel could be correlated due to spatial impact. Specially, 
in the sectorized celluar network outdoors, signals from the 
mobile users in a sector could be distributed within a narrow 
angle ]. The Vandermonde matrix can be used to depict the 
channel model of the outdoor LOS propagation path for the 
large antennas system 000 . Most of the current 
research on massive-MIMO is limited to theory based on the 
analytical and simulation results, and lacks of the experimental 
support. Thus our experimental research is motivated to verify 
the theoretical results and to find the effect of massive MIMO 
in real propagation environment. We built a prototype of 



Fig. 1: System model of MU-MIMO system with TVt-antenna 
base station and N r mobile users 


massive MIMO of testbed with the commercial SDR platfrom 
- USRP N210. With 30 high performance computer and 
70 SDR nodes, we developed a massive MIMO experiment 
testbed supporting up to 30 x 30 scheme. We also developed 
the software to support the distributed data collection and 
computing architecture. Our preliminary experimental work 
based on this testbed include 

• Implement of the TDMA protocol to ensure the orthog¬ 
onality in channel measurement; 

• Compute the massive MIMO capacity based on the 
measured channel matrix, pre-coding is also considered; 

Our experimental results show the linearity of the capacity 
with the increasing number of antennas by assuming it is the 
same with the number of mobile users. Further, more impor¬ 
tantly, we found out that our measured capacity agrees with 
the simulated capacity under Vandermonde channel model. It 
is the first time that the feasibility of Vandermonde channel 
model is demonstrated by experimental result. 

tntech 

II. Theoretical Model 

Considering a point-to-point MIMO system in which there 
are N t transmit antennas and N r receive antennas. For the 
narrow band channel, assuming the transmitters transmit the 
data simultaneously, the mathematical description of the signal 
model is as below: 
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y = v/pHs + n (1) 

where s G C Nt is the transmitted vector satisfying the 
E [ss H ] = 1, p is the common SNR(signal-to-noise-ratio) 
at each receive antenna and n is a vector of complex 
Guassian noise complying to zero mean and unit-variance 
i.i.d(independent and identical distribution), H G £ N r* N t is 
the channel response matrix. 

A. Capacity and achievable rates 

Under the assumption that the receiver has perfect knowl¬ 
edge of channel matrix H, the capacity of the N t x N r MIMO 
channel is computed by 

C = log 2 det (r Vr + £ HHff ) (2) 

where Uv r is the identity matrix and the H means hermitian 
transposition. 
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Fig. 2: The experiment deployment of the MU-MIMO system. 


B. Vandermonde Channel Model 

In (5) (6) (7) the Vandermonde random matrix 
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e -j27r^ sin(0jv) 
e -j27r(M-l)f sin (0 N ) 

(3) 

is introduced to describe the channel model for a base station 
receiver with M antennas and N mobiles, where d is the 
antenna spacing and A is the wavelength. As in fig. |TJ the 
angles of arrival are supposed to be uniformed distributed 
within (—a, a). The elements of the Vandermonde matrix can 
also have phases with uniform distribution for comparison. 
The received signal at the base station is given by 

y = VP^s + n (4) 

where y, s, n are respectively the Mxl received vector, the 
N x 1 transmit vector, and the Mxl additive noise, V is 
the Vandermonde channel matrix, P is the power gain matrix 
which can be set as identity matrix in simulation. 

The Vandermonde model, compared with Gaussian channel 
model, is close to the real massive MIMO system from the 
architecture point of view. 
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III. System design and experiment setup 


A. Channel Measurement 


The channel matrix is key to compute the MIMO system 
capacity. Denote hij as the frequency response for the link be¬ 
tween the transmit antenna i and receive antenna j. Assuming 
the narrowband channel, the channel can be measured by the 
PN sequence correlation u 

At the transmitter, a PN sequence with length of 512, 
x = [cq, • • • Q, • • • C 510 ], is transmitted. The auto-correlation 
function of the PN sequence is 


Rxx (jl) — 


u 


n = 0, TV, —IV, 2AT, -21V, 
otherwise 


(5) 


At the receiver, the correlation between the received se¬ 
quence, y and the known transmitted PN sequence, x is 

Rxy (jt) = hijRxx (v) ( 6 ) 

Thus the h tJ is computed based on eq. [6] 

In the real channel measurement, the TDD(Time division 
duplex) is introduced to measure the channel matrix by ensur¬ 
ing the orthogonality of the channel sounding signal. 

B. Experiment Deployment and System Architecture 

In the wireless networking systems lab (WNSL) of Ten¬ 
nessee Tech university (TTU), we built a massive MIMO 
testbed based the general SDR platform - USRP. Seventy 
USRPs controlled by thirty high performance PCs are de¬ 
ployed on the 4th floor of Clement Hall on the campus of 
TTU. As shown in Fig. [2| two groups of USRPs are placed 
in two separate rooms. Forty USRPs are deployed in the 
WNSL to emulate multi-antenna base station, while the rest 
of the USRPs are randomly arranged in another room acting 
as mobile users with single antenna. 

The massive MIMO system is composed of data plane 
and control plane, as depicted by Fig. [3] The data plane 
is responsible for the channel measurement and waveform 
data processing. On the data plane, 30 PCs are configured 
as baseband processing units. Every baseband processing unit 
is connected to arbitrary number of USRP front-ends acting 
as either transmitter or receiver, through Ethernet. On the 
control plane, one central control PC is connected with all 
the baseband processing units to control the procedures of 
the experiment including the start and stop of measurement, 
parameters configuration, etc. The communications between 
the baseband processing unit and control unit is also over 
Ethernet and the socket/TCP programming. 

C. Massive MIMO testbed over time division protocol 

The channel matrix measurement is a challenging task that 
is achieved by our massive MIMO testbed. Firstly, we want to 
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Nr Receivers 



over USRP SDR. 


measure the data within the time frame that the channel can 
be regarded as static. Secondly, there are so many links (30 x 
30) to be measured within this time frame. The orthogonality 
needs to be ensured that the measured result for each link is 
not interfered by all the rest of the links. We choose the time 
division method to separate the channel measurements for all 
the links. Under the time division strategy, our massive MIMO 
testbed is equivalent to a TDMA network. 30 transmitters only 
send the PN sequence within their own time slots, while the 
mobile users can be set as receive only mode. 

During the experiment, suppose the number of the trans¬ 
mitters is the same with the number of mobile users, so we 
actually deployed a 30 x 30 MIMO system which is regarded 
as a massive MIMO testbed as the scale is much bigger than 
the current commercial MIMO system. 

D. Synchronization 

The timing synchronization for all the TDMA nodes is 
achieved by the GPS signal as shown in Fig. [3] The timing 
precision provided by the GPS signal reaches nano-second, 
which is sufficient for the time slot scheduling. Every USRP 
is equipped with a dedicated GPS antenna that can feed the 
timing information to the front end. The baseband samples 
sent from the PC to the USRP can be arranged within the 
specified time frame for transmitting by adding the time tag 
information at the start of sample series. 

We also provide the clock synchronization for all the USRP 
front-ends by the clock distributing board AD9523 as in 
figure. [3] The common clock between the transmitters and 
receivers helps to remove the carrier frequent offset(CFO) that 
is introduced by different clock sources. 


system. 30 x 30 configuration costs 1.55 to get a round of data. 
Basically the channel can be regarded as static within several 
seconds for our current static deployment because there no 
other impact during the experiment. 

We did the simulations for the capacity under the Vander¬ 
monde channel model, by adjusting two parameters d and a 
based on section II-B Fig. [4] shows the experimental results 
compared with the simulations. The experimental frequency 
is 926MHz, the corresponding wavelength is A = 32.4cm. 
The average distance between the adjacent nodes for base 
station is about 20cm. We selected the parameter d in 
{1.0A/2,1.1A/2,1.2A/2,1.25A/2} and did the simulation for 
Vandermonde model by adjusting a for each fixed d. 

Firstly, we demonstrate the linearity of the capacity with 
the number of antennas by the measured channel capacity in 
Fig. [4] The capacity of the system increases linearly with the 
number of antennas, for both the simulated Gaussian channels, 
Vandermonde channels and the measured channels. 

Secondly, we found the feasibility to fit the massive MIMO 
channel model to the Vandermonde channel model, by the 
comparison of capacities between the real channel measure¬ 
ment and the Vandermonde channel model. Fig. [4] shows 
observations including : 

• The capacity of Vandermonde channel decreases when a 
is getting smaller for fixed d. 

• The capacity of Vandermonde channel increases when d 
is bigger for fixed a. 

• For every fixed d , we can find an a to fit the Vandermonde 
model to measured channel in capacity with minimum 
mean error (MME). Table. [I] shows the MME value for 
different a under each fixed d. Bold a represents the 
optimal a for each d. 

We get insight from above that the Vandermonde channel 
model is more realistic than the Gaussian channel model in 
depicting the massive MIMO system. The Gaussian model 
does not reflect the topology of the real MIMO system while 
the Vandermonde model is more illustrative as the parameter a 
and d are contained. Also, the Gaussian model always has the 
ideal capacity which is hard to be achieved in the real system. 
We can always find an equivalent Vandermonde model for the 
measured massive MIMO channel by our real experimental 
data, when their capacities are very close under optimal pair 
of a and d. 

In Table. [TJ the parameters pair (d, a) with values of 
(1.25A/2, 28 degree) gives the minimum MME value. This 
optimal parameters pair are very close to that in the real 
deployment (20cm, 26 degree). 


a 

MME-di 

a 

MME-^2 

a 

MME-^3 

a 

MME - d 4 

35 

40.7338 

35 

86.2465 

30 

52.0346 

30 

77.9740 

30 

98.9845 

32 

39.7652 

29 

39.7046 

28 

37.9838 

26 

189.545 

30 

51.8876 

26 

66.8592 

26 

50.9662 


IV. Experiment results 

This massive MIMO testbed can be scaled to at most 30 x 
30. The measured channel matrix is used to calculate the real 
MIMO system capacity. 

In the experiment, if we set the time slot duration to 50 ms, 
we can get a round of TDMA data by 200m5 for 4 MIMO 


TABLE I: Minimum mean error for different pair of d and a. 

V. CONCLUSION 

This paper presents the design and implementation of a 
novel SDR based massive MIMO testbed using TDMA proto- 
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d = 1.0’lambda/2 



d=1.1*lambda/2 



(a) d = 1.0A/2 


d=1.2*lambda/2 



(b) d= 1.1A/2 


d=1.25*lambda/2 



(c) d= 1.2A/2 


(d) d= 1.25A/2 


Fig. 4: Comparison of measured channel with Vandermonde channel model by adjusting d and a , in terms of capacity. 


col. The MIMO channel matrix is measured through the exper¬ 
iment with the prototype of the testbed. The system capacity 
is computed from the measured channel matrix. The testbed 
is validated by the linearity of capacity with the increasing 
number of antennas. We also found that the measured channel 
capacity agrees with that of the Vandermonde channel model 
with the optimal parameters: the range of the arrival angles and 
the base station antenna distance which are very close to those 
in experimental deployment. It is the first time the feasibility 
of the Vandermonde channel model is shown in the real 
experimental result. Based on our work, we recommend the 
Vandermonde channel model to be used in future theorectical 
research, as it is more realistic than the current widely used 
Gaussian model. 
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